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Based on recent developments by the authors a next-to-leading order spin(l)-spin(2) Hamiltonian 
is derived for the first time. The result is obtained within the canonical formalism of Arnowitt, 
Deser, and Misner (ADM) utilizing their generalized isotropic coordinates. A comparison with 
other methods is given. 
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I. INTRODUCTION 

In this paper we present the next-to-leading (NLO) order gravitational spin(l)-spin(2) dynamics in Hamiltonian 
form. The result is based on the ADM canonical formalism [l[ for spinning classical objects recently derived by 
Steinhoff, Schafer, and Hergt which already has shown its power by the derivation of the Hamiltonian of two 
spinning compact bodies with next-to-leading order gravitational spin-orbit coupling, lately obtained by Damour, 
Jaranowski, and Schafer [3|. 

The following notations will be used throughout the paper: The canonical spin tensor of the a-th object (particle) 
is the euclidean spin tensor S l J = S a ij = e^ k S k , S a = [S k ) = (S a k), and it also holds by definition S a ij — e-ik^jiSaki 
with eik being the symmetric root of the symmetric 3-metric 7^. S a %j and S a ij fulfill the conserved-length relations 
SaijJ lk '~f^ 1 S a ki = S a ij Saij = const. The mass parameter of the a-th particle is denoted m a . The 4- vector is 
the unit vector orthogonal to the spacelike hypersurfaces t = const.; its components are — (-N, 0,0,0), where 
N is the lapse function. The short-cut notation —np a is used for yj m a + j^PaiPaj , where (p a i) = p a denotes the 
canonical momentum of the a-th particle. The canonical particle position variables are (x a ) — x a and the velocities 
read v a = (v a ) — {x l a ), where the dot means coordinate time derivative. 

Our units are c = 1, where c is the velocity of light. G will denote the Newtonian gravitational constant. Greek 
indices will run over 0, 1, 2, 3, Latin over 1, 2 if from the beginning of the alphabet and 1, 2, 3 if from the middle. For 
the signature of spacetime we choose +2. 



II. SPINNING OBJECTS IN THE ADM FORMALISM 
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Recently in it has been shown that the matter source parts of the energy and momentum constraint equations, 
respectively 7i M and , are given in terms of canonical position, momentum, and spin variables in the form, to the 
post-Newtonian orders indicated, 
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These expressions are sufficient for the derivation of the Hamiltonian in [3J] . The applied equal-time Poisson brackets 
read 



{xi,p aj } = % , {S l a j ,S kl } = S lk S a l - S Jh SS - 5 a Si k + S^Sf , zero otherwise . 



(2.3) 



It is important to note that the expressions for Ti. M and Tif 1 are given in the ADM transverse-traceless (ADMTT) 
gauge which refers to generalized isotropic coordinates defined by the conditions 



Hj = ( 1 + \<t> ) 5ij + h 
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with hj^ being transverse and traceless (hj^j = 0, hj^ = 0). The square root of the metric takes the form, to 
sufficient approximation for the Hamiltonian, 

e« = U + \^j Sii+O^h^) . (2.5) 

In this approximation, Sj = S ai a l3 = S a ijj h = S l aj = S alJ = S% holds. 
The canonical conjugate to hj^ , lg * G 7r^ T , i.e., 

{^ T (x,i),7r*; i T (x , ,t)} =^ Tfez (x-x'), zero otherwise , (2.6) 



16ttG 



will play no role in the calculations of the present paper. Only the longitudinal part of , tt 1 ^ , contributes (notice: 
w ij _ ^ij _|_ tt^j-.). Thus, the ADM Hamiltonian H will not depend on 7r^! T , i.e., 

H[x a ,p a ,hJ?} = J d 3 .TA0. (2.7) 

Therefore, hj? is allowed to be replaced by matter variables (after solution of the evolution equations which in 
our approximation is an elliptic equation only) to get an autonomous Hamiltonian. Otherwise the transition to a 
Routhian would have to be performed, see jij. The equations to be solved for the obtention of the autonomous matter 
Hamiltonian read, on the one side 



7iklji7r ij TT kl + -7y7«^'7r W + V^R(3) = 1GkGH m , (2.8) 



-li j -K jk . k = ^GHf , (2.9) 

where i?( 3 ) denotes the Ricci scalar of the t = const, slices, 7 is the determinant of the 3-metric, and ; the 3-dim. 
covariant derivative, and on the other side 

= {H,tt 1 4 t } + O(1/c 6 ). (2.10) 
After tedious calculations, the NLO spin(l)-spin(2) interaction part of the Hamiltonian results in 

H SS° = o ^— 3-[§((Pi x S i) ' n i2)((P2 x S 2 ) • n 12 ) + 6((p 2 x Si) ■ n 12 )(( Pl x S 2 ) • n 12 ) 

- 15(Si • ni 2 )(S 2 ■ n i2 )(pi • ni 2 )(p 2 • n i2 ) - 3(Si • ni 2 )(S 2 ■ n i2 )(pi • p 2 ) 

+ 3(Si • p 2 )(S 2 • n i2 )(pi • n i2 ) + 3(S 2 • pi)(Si • ni 2 )(p 2 ■ n i2 ) + 3(Si • pi)(S 2 • n i2 )(p 2 • n i2 ) 
+ 3(S 2 ■ Pa )(Si • ni 2 )(pi ■ n 12 ) - |(Si ■ p 2 )(S 2 • Pl ) + (Si ■ Pi)(S 2 ■ p 2 ) 

- 3(Si • S 2 )(pi • ni 2 )(p 2 • ni 2 ) + |(Si ■ S 2 )(pi • p 2 )] (2.11) 
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+ a 3 [~((pi x Si) • n 12 )(( Pl x S 2 ) • n 12 ) + (Si • S 2 )(pi • n 12 ) 2 - (S x • ni 2 )(S 2 • Pi)(pi ■ n 12 )] 
2m 1 r 12 

3 

[-((P2 x S 2 ) • ni 2 )((p 2 x Si) • H12) + (Si • S 2 )(p 2 • m 2 ) 2 - (S 2 • ni 2 )(Si • p 2 )(p 2 • n i2 )] 



2m\r\ 2 



+ 6(mi + m2) [(S! • S 2 ) - 2(Si ■ ni 2 )(S 2 • m 2 )] , 

where ri 2 = |xi — x 2 | is the euclidean distance between the two particles and ni 2 denotes the unit vector ri 2 ni 2 = 
xi - x 2 . 

III. DIFFERENT DERIVATION OF THE SPIN(1)-SPIN(2) HAMILTONIAN 

We follow here the procedure described in Q. The implementation of spin into the Eq. (4.9) of [H results in 

ispi„ = _ y ( 3mAjj 2S \ Kb (3 1} 
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The NLO spin(l)-spin(2) interation Hamiltonian is given by 



H, 



NLO _ o ci nj _ o s P in ( 2 ) « — n 8 ?"^ 1 ) « 

SS — "(4)ij S l ^2 ~ 5J (4) ' 3 1 — 1J (4) ' 3 2 
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Using Eq. (4.10b) in Ref. @], but calculated for metric functions resulting from our matter source terms, the obtained 
Hamiltonian coincides with the one calculated in the present paper. 

In this completely independent approach also the matter source part h^l 7 Tkl of the evolution equations con- 
tributes, 
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and lapse and shift functions have to be determined too. Details of the calculations can be found in (2j. 



IV. CONSISTENCY 



The correctness of the derived spin(l)-spin(2) Hamiltonian can best be verified by the construction of the global 
Poincare algebra. The generators of the global Poincare algebra are the total linear momentum P = Pi + P2, the 
total angular momentum J = Xi x pi + x 2 x p2 + Si + S2, the total Hamiltonian H, and the total center-of-mass 
generator G. The latter object is defined by G = — 16 * G J d 3 xxA0, e.g., see and turns out to be 



G = G 
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(S 2 ■ n 12 )Si - (Si • m 2 )S 2 + (3(S 2 ■ n 12 )(Si ■ n 12 ) - (Si ■ S 2 )) 



xi + x 2 

ri2 



(4.1) 



where Gso denotes the NLO spin-orbit coupling contribution as given in Q and Gpm is the point-mass part from 
6] . It is straightforward to show that the above generators do fulfill the Poincare algebra. 



COMPARISON WITH OTHER METHODS AND RESULTS 



In a recent paper based on Ref. Q, Porto and Rothstein derived a next-to-leading order spin(l)-spin(2) potential 
using an action approach with spin supplementary conditions imposed on the action level. A consistency calculation, 
however, which would have shown that the spin supplementary conditions are preserved under the variational principle, 
and thus under the equations of motion, has not been undertaken nor has it been shown that, in contrast to their 
claim, the used position, velocity, and spin variables are those that relate to canonical ones in standard manner, cf., 
e.g., [9(. A consistency check of their intuitive canonical approach by an independent method is therefore necessary, 
which is the subject of this Section. 

The relevant part of the Lagrangian of Porto and Rothstein reads I/ PR = ^mivf + \"m,i^f\ — Vn — Vgo — ^ss R wnere 
Vn is the Newtonian potential, Vg S is given by Eq. (12) in and the leading order spin-orbit coupling potential 
function VgQ results from, e.g., [3[, 



v so — -1 



3 3 

-m 2 ((vi x Si) • m 2 ) - 2m 2 ((v 2 x Si) • n X2 ) - -mi((v 2 x S 2 ) • n i2 ) + 2m x ((vi x S 2 ) • n i2 ) 



(5.1) 



The canonical momenta are given by p = — , i.e., 



Pi = mivi 
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p 2 = m 2 v 2 + — r 

r 12 



-m 2 (Si x n i2 ) + 2mi(S 2 x n i2 ) 



-mi(S 2 x m 2 ) + 2to 2 (Si x n 12 ) 
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(5.2) 
(5.3) 



The Hamiltonian of Porto and Rothstein then takes the form H — Hn + -ffgo + Hg§. We can get H by replacing 
the velocities by canonical momenta in the following expression 



ff PR = vi. P i+v 2 .p 2 -L PR . 



(5.4) 
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Note that the -g^ 2 - terms do not contribute. The difference between -ffgg LO and reads 



G 



-[3(Si • P2 )(S 2 • ni 2 )(pi • ma) + 3(S 2 ■ Pl )(Si ■ n 12 )(p 2 • n 12 ) 



2miTO 2 rJ 2 1 

- 2(Si • P2 )(S 2 ■ pi) - 6(Si ■ S 2 )(pi • ni 2 )(p 2 ■ n 12 ) + 2(Si ■ S 2 )(pi • p 2 ) 
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[3(Si • S 2 )(pi ■ n 12 ) 2 - (S : ■ S 2 )pf - 3(Sa • n 12 )(S 2 • Pl )( Pl • n 12 ) + (S 1 ■ Pl )(S 2 • Pl )] 
[3(Si • S 2 )(p 2 • n 12 ) 2 - (Si • S 2 )p 2 - 3(S 2 • n 12 )(S 1 • p 2 )(p 2 • n 12 ) + (Si • p 2 )(S 2 • p 2 )] 



(5.5) 



Obviously, there is agreement at O (G 2 ). 

There should exist an infinitesimal generator g for a canonical transformation such that SH^g° = {Hn, g}. Plugging 
in the ansatz 
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and comparing O (G) terms gives 



The vanishing of O (G 2 ) terms yields 



a = 0. 



TO2 



c= -, d = 0. 
2 



-(Si • ni 2 )(S 2 • ni 2 )(p 2 • ni 2 ) 



-a + b-d = 0. 



(5.7) 



(5.8) 



which is incompatible with the canonical transformation that is needed for the O (G) terms. 

In Ref. which is a short reply to the first version of the present paper [TTj] . Porto and Rothstein pointed out 
that their result in 0] is incomplete in the sense that it only includes contributions from spin-spin diagrams. It was 
not realized in 0] that spin-orbit diagrams also contribute to the next-to-leading order spin(l)-spin(2) interaction. 
If these contributions are included, the canonical transformation defined by (|5.6p and (|5.7p leads to an additional 
agreement at O (G 2 ), and thus to full agreement. It should be recalled that our derivation is completely different and 
includes all contributions to the spin(l)-spin(2) interaction from the very beginning. 



Acknowledgments 



This work is supported by the Deutsche Forschungsgemeinschaft (DFG) through SFB/TR7 "Gravitational Wave 
Astronomy" . 



[1] R. Arnowitt, S. Deser, and C. W. Misner, in Gravitation: An Introduction to Current Research, edited by L. Witten (Wiley, 

New York, 1962), p. 227, |arXiv:gr-qc/0405109| 
[2] J. Steinhoff, G. Schafer, and S. Hergt, The ADM canonical formalism for gravitating spinning objects, Phys. Rev. D, 

submitted. 

[3] T. Damour, P. Jaranowski, and G. Schafer, Phys. Rev. D 77, 064032 (2008). 
[4] P. Jaranowski and G. Schafer, Phys. Rev. D 55, 4712 (1997); 57, 7274 (1998). 
[5] T. Regge and C. Teitelboim, Ann. Phys. (N.Y.) 88, 286 (1974). 

[6] T. Damour, P. Jaranowski, and G. Schafer, Phys. Rev. D 62, 021501 (R) (2000); 63, 029903(E) (2001). 
[7] R. A. Porto and I. Z. Rothstein, Phys. Rev. Lett. 97, 021101 (2006), |arXiv:gr-qc/0604099fr rl. 



5 



[8] R. A. Porto, Phys. Rev. D 73, 104031 (2006). 
[9] K. Yee and M. Bander, Phys. Rev. D 48, 2797 (1993). 
[10] R. A. Porto and I. Z. Rothstein, Comment on 'On the next-to-leading order gravitational spin(l)-spin(2) dynamics' by J. 

Steinhoff et al. larXiv:0712.2032V l [gr-qc]. 
[11] J. Steinhoff, S. Hergt, and G. Schafer, On the next-to-leading order gravitational spin(l)-spin(2) dynamics, 
larXiv:0712.17 16vl [gr-qc]. 



